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We report results on the Coulomb-gauge ghost propagator and the color-Coulomb potential com- 
puted in two lattice gauge-field ensembles: (1) configurations derived from our recently proposed 
Yang-Mills vacuum wave functional in 2 + 1 dimensions, and (2) lattices generated by Monte 
Carlo simulations of the three-dimensional Euclidean SU(2) lattice gauge theory with the Wilson 
action. We observe remarkable agreement between the ghost propagators in both ensembles, but 
some differences in the potentials. Those originate from rare configurations with very small val- 
ues of the lowest eigenvalue of the Coulomb-gauge Faddeev-Popov operator. If the same cuts on 
such exceptional configurations are applied in both ensembles, then the color-Coulomb potentials 
are also in reasonably good agreement. 
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1. Proposal for an approximate vacuum wave functional 



Confinement is supposed to be encoded in properties of the vacuum of quantized non-abelian 
gauge theories. In the hamiltonian formulation in D = d + 1 dimensions and temporal gauge, the 
vacuum wave functional satisfies the Schrodinger equation: 
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together with the Gauss-law constraint: 
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At large distance scales one expects that the wave functional assumes the effective form: 
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It has a property of dimensional reduction [|], ||, |3j] : The computation of a spacelike loop in 
d + 1 dimensions reduces to the calculation of a Wilson loop in Yang-Mills theory in d Euclidean 
dimensions. However, the true vacuum wave functional cannot be just like this - it implies not 
only Wilson's area law, but also e.g. exact Casimir scaling of higher-representation string tensions, 
which is not observed at asymptotic distances. 

Recently, we have proposed a simple approximation to the vacuum wave functional of the 
(2+ 1) -dimensional SU(2) Yang-Mills theory jH] 1 : 
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where B a {x) = F" 2 {x) denotes the color magnetic field strength, ^[A] is the covariant derivative in 
the adjoint representation, 3 2 = 2$k ■ S'k the covariant laplacian in the adjoint representation, Xq is 
the lowest eigenvalue of (— 3) ), and m is a constant (mass) parameter proportional to g 2 ~ 1//3. 



The expression (1.4) is written in the continuum notation, but assumed to be properly defined on a 
lattice, where we choose 
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and (x) denotes the link matrix in the fundamental representation. 



(1.5) 
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2. Arguments in favor of the proposed vacuum wave functional 

In the original paper we have supplied a series of arguments supporting the proposed form 
of the Yang-Mills vacuum wave functional: 



'A similar form, but without the - crucial in our opinion - subtraction of Xq, was proposed by Samuel \a\. 
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1. In the free-field limit (g — > 0), *Po[^] becomes the well-known vacuum wave functional of 
electrodynamics. 

2. The proposed form is a good approximation to the true vacuum also for strong fields constant 
in space and varying only in time. 

3. If we divide the magnetic field strength B(x) into "fast" and "slow" components, the part of 
the vacuum wave functional that depends on B s i ow takes on the dimensional-reduction form. 
The fundamental string tension, at a given j8, is then easily computed as Op = 3m/4j8. 

4. If one takes the mass m in the wave functional as a free variational parameter and computes 
(approximately) the expectation value of the Yang-Mills hamiltonian, one finds that a non- 
zero (finite) value of m is energetically preferred. 



3. Lattice evidence 



The above (analytic) hints are encouraging, but in no way sufficient to convince anybody that 



the simple vacuum wave functional, Eq. (1.4), is close to that of the true Yang-Mills ground state. 
To assess how good or bad the approximate state is, some numerical tests are inevitable. For that 
purpose, we compared a set of quantities computed in two ensembles of lattice configurations: 

1 . "Recursion " lattices — independent two-dimensional lattice configurations generated with 
the probability distribution given by the proposed vacuum wave functional, with m fixed at 
given j6 to get the correct value of the fundamental string tension Of(J3)- The recursion 
method was described with all details in Ref. [Q]. 

2. Monte Carlo lattices — two-dimensional slices of configurations generated by Monte Carlo 
simulations of the three-dimensional euclidean SU(2) lattice gauge theory with the standard 
Wilson action; from each configuration, only one (random) slice at fixed euclidean time was 
taken. 

3.1 Mass gap 

The first such test was performed in Ref. [Q]. We computed the equal-time connected B 2 -B 2 
correlator and determined the value of the mass gap from a best fit to its exponential fall-off at large 
distances. The result for recursion lattices is compared in Fig. [I] to the values of the + glueball 
mass computed in high-statistics simulations of the three-dimensional Yang-Mills theory by Meyer 
and Teper [Q]. The deviations are at the level of at most 6%. 

3.2 Coulomb-gauge quantities 

Another set of quantities of interest are defined in the Coulomb gauge. It was argued by Gri- 
bov [f7|] and Zwanziger [Jsj] that the low-lying spectrum of the Faddeev-Popov operator in Coulomb 
gauge probes properties of non-abelian gauge fields that are crucial for the confinement mechanism. 
The ghost propagator in Coulomb gauge and the color-Coulomb potential are directly related to the 
inverse of the Faddeev-Popov operator, and play a role in various confinement scenarios. In par- 
ticular, the color-Coulomb potential represents an upper bound on the physical potential between 
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Figure 1: Mass gaps from ensembles of recursion lattices at various lattice couplings, compared to the + 
glueball masses in D = 2 + 1 obtained in Ref. [||] via standard lattice Monte Carlo methods. (From Ref. [Q].) 

a static quark and antiquark, which means that a confining color-Coulomb potential is a necessary 
condition to have a confining static quark potential [B]. 



Our aim [ |10| ] was to see how well the proposed vacuum wave functional can reproduce the 
values of Coulomb-gauge observables that can be obtained by standard lattice MC techniques. For 
that purpose we had to measure them in Coulomb-gauge configurations distributed with the prob- 
ability following from the square of our temporal-gauge vacuum wave functional. In the operator 
formalism, the minimal Coulomb gauge is a gauge fixing within the temporal gauge of the remnant 



local gauge invariance. (See Ref. Qll| ] for a detailed discussion of this point.) The wave-functional 
in Coulomb gauge is the restriction of the wave functional in temporal gauge to transverse fields in 
the fundamental modular region A: 

^Coulomb [ A j = ^ A A± g A (31) 

The vacuum expectation value of an operator Q in Coulomb gauge can then be computed from 

(Q) = (^"lelA^]!^ ") = OPolG [ a A] pF >, (3.2) 

i.e. we generate configurations following the probability distribution SPq, transform them to the 
Coulomb gauge, and evaluate the observable Q in the transformed configuration. (£2 denotes the 
gauge transformation that brings the configuration A to the minimal Coulomb gauge.) From the 



path-integral representation of the vacuum state, we may also go from ( p.2j ) to 



(Q) = (Q 



t ) 



(3.3) 



where the right hand side is the expectation value obtained in D = 3 Euclidean dimensions, and CI 1 
is the gauge transformation which takes the gauge field ona( = to time-slice into Coulomb gauge. 

Fig. ^displays results for recursion and Monte Carlo lattices at j8 = 9 on 32 2 lattice. The ghost 
propagator in Coulomb gauge was computed from the inverse of the Faddeev-Popov operator (in 
the subspace orthogonal to trivial constant zero modes due to lattice periodicity) 
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Figure 2: The Coulomb-gauge ghost propagator (a) and the color-Coulomb potential (b) at j3 = 9 on 32 2 
lattice. 
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with ^(jc) = ^Tr[J7i:(^)] and a L k (x) = jjTr[o c Uk(x)], while the color-Coulomb potential between a 
static quark and antiquark located at points x and y is proportional to 



V(R) = -/(^[A]-\-V 2 )^[A]- 1 ) 
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The agreement of the ghost propagator computed in both sets of lattices is almost perfect 
[Fig. ||(a)], the differences are at the level of the size of symbols. On the other hand, there is a 
considerable deviation of color-Coulomb potentials computed for the MC ensemble from those for 
recursion lattices [Fig. |2|(b)]. Fortunately, it is not difficult to explain the origin of this deviation. 
There exist, in both ensembles, "exceptional" configurations with a very low (though still positive) 
lowest nontrivial eigenvalue of the Faddeev-Popov operator. These configurations were extremely 
difficult to gauge-fix to the Coulomb gauge. If one evaluates the potential in each single configura- 
tion, the exceptional ones possess a very high absolute value of the potential at the origin, |V(0)|. 
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Table 1: Numbers of configurations in different ranges of |V(0)| at J3 = 9 on 32 2 lattice. 
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Figure 3: The color-Coulomb potential for K = 10 (left) and K = 20 (right). J3 = 9, 32 2 lattice. 
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Figure 4: The color-Coulomb potential for K = 100 (left) and without cut (right). j3 = 6, 24 2 lattice. 



One can then classify configurations by their values of |V(0)|, and evaluate average potentials from 
sets of configurations satisfying a number of cuts: {|V(0)| < fc ; -, i = 1,2,... ,K}. For illustration, 
Table [I] lists numbers of configurations in a couple of | V (0) | bins at j8 = 9 (32 2 lattice). 

Fig. U shows results for two values of fc at j3 = 9 (32 2 lattice). The potentials agree quite 
well for lower fc cut (satisfied by about 80% lattices), but the differences (and errorbars) grow with 
increasing K. The problem appears even more spectacular at j3 = 6 (24 2 lattice, Fig. There 
is close agreement of color-Coulomb potentials between Monte Carlo and recursion lattices up to 
the cutoff k as high as 100, but then a single recursion lattice with extremely high value of |V(0)| 
completely distorts the picture. Generally, the intervals with high |V(0)| values are very unequally 
populated and lead to wild fluctuations of the results, when the cut K is increased. We believe 
that approximate agreement would be restored with sufficient (but obviously huge) statistics, even 
though some differences might persist. 



4. Conclusion and outlook 



The proposed vacuum wave functional for the temporal-gauge SU(2) Yang-Mills theory in 
2+1 dimensions, Eq. (1.4), seems a fairly good approximation to the true ground state of the theory. 
We have added two new pieces of evidence in its favor: 
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1. The ghost propagator in Coulomb gauge is practically identical in recursion and Monte Carlo 
ensembles. 

2. With the same statistics of rare exceptional configurations we expect also the color-Coulomb 
potential from recursion lattices to be close to that determined from Monte Carlo lattices. 

Our further goals are to compare consequences of our proposal to others existing in the liter- 



ature (see e.g. Q12|]), to determine the wave functional in numerical simulations for typical field 



configurations by the method of Ref. [13], to improve on the variational estimate of the parameter 



m, and to generalize the proposal to the realistic case of 3 + 1 dimensions. 
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